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Abstract 

In the paper we define the convergence of compact fuzzy sets as a convergence of ct-cuts in 
the topology of compact subsets of a metric space. Furthermore we define typical convergences of 
fuzzy variables and show relations with convergence of their fuzzy distributions. In this context 
we prove a general formulation of the Strong Law of Large Numbers for fuzzy sets and fuzzy 
variables with Archimedean i-norms. Next we dispute a structure of fuzzy logics and postulate 
a new definition of necessity measures. Finally, we prove fuzzy version of the Glivenko-Cantelli 
theorem and use it for a construction of a complete fuzzy measurement theory. 

1 Introduction 

1.1 Motivations 

The application of fuzzy set theory in the description of empirical data becomes more popular in the 
recent years. A consistent theory requires algorithms which allow to construct a membership function 
of a given process with use of empirical data. This problem is completely solved in the theory of 
probability where by the Strong Law of Large Numbers a probability has an interpretation of frequency 
and a convergence of process estimators is guaranteed by Central Limit Theorem and Glivenko- 
Cantelli theorem. These theorems constitute a basis for further statistical analysis, construction of 
estimators and hypothesis verification. 

The idea of fuzzy variables, counterparts of probabilistic random variables, reaches papers of 
Nahmias |17j and Rao and Rashed [20] . However, analogous theorems for fuzzy sets or fuzzy variables 
are not general enough. The most important and applicable result on the Law of Large Numbers is 
due to Triesch [23) . However, this result was never interpreted in the language of fuzzy variables. 

Before laws of large numbers may be proved, different types of convergence of fuzzy sets and fuzzy 
variables must be discussed. This requires to define at least a topological space of all fuzzy sets and 
fuzzy variables. Some general results on convergence of fuzzy sets are due to Diamond and Kloeden 
[4] and Kaleva [12]. A convergence of fuzzy variables was used implicitely in the papers of Fuller 
[B], Triesch [S3] and Williamson 27J, but a general definition and dependencies between different 
convergences were not studied in detail. 

Laws of Large Numbers were introduced in a few different versions. In the papers of Hong [10j . 
Hong and Ro [TT] and Triesch [52] a pointwise convergence of fuzzy intervals is discussed. These 
results are based on analytical expressions for sums of fuzzy intervals developed by Fuller [7] , Markova 
[T5] . Mesiar [IB], Hong [pj, Hong and Hwang [5] and other which are valid only in special cases for 
example when their membership functions of the slopes are logarithmically concave. On the other 
hand a method of building up a membership function from a sample based on probability-possibility 
transformations discussed in the paper of Dubois et al. [3J frequently leads to convex membership 
functions. 
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A different approach to the Law of Large Numbers is presented in the papers of Fuller [6] and 
Triesch [23]. We will show that a convergence analyzed in these papers is in fact a convergence of 
fuzzy variables. We will give an alternative proof of the theorem 3 in [2*5] . 

A discussion on convergence leads to a fuzzy counterpart of Glivenko-Cantelli theorem which we 
will proof in the last section. This theorem, crucial for a fuzzy data analysis, was never proved before. 

In this paper we construct a complete theory of fuzzy variables. We discuss convergence of fuzzy 
variables and give a proof of the Strong Law of Large Numbers. Next, we construct an estimated 
membership function of a given process from a series of empirical data similarly to the procedure 
well-known from probability theory. Our fuzzy sets represent the whole process rather than a simple 
measurement result and describes both components of uncertainty, a random and a systematic one. 
This procedure does not require probability-possibility transformations or any probabilistic interpre- 
tation of fuzzy data. 

Our philosophy is to look at fuzzy sets as a collection of a-cuts rather than a membership function. 
Hence we would like to study a convergence of fuzzy sets as a convergence of their a-cuts. This requires 
to define a topology on the space of a-cuts which can be naturally done in case of compact fuzzy sets. 
In this fashion we want to estimate a-cuts directly from empirical data rather than a membership 
function of a given fuzzy process and show a convergence of such a procedure to a fuzzy interval 
representing this process. In such an approach we use positional statistics such as modal value and 
median which are more natural in this context. 

Finally, we must underline that our fuzzy variables are not fuzzy random variables. The results 
for fuzzy random variables are similar to those known from probability and were intensively studied. 
Laws of Large Numbers are well-known in this case, e.g. [2*1] . [13] . 

1.2 Definitions 

We use the following notations. By R we denote real numbers, N = {0, 1,2,.. .}, |~— ] and [— J stand 
for floor and ceil functions respectively. For any set ft by 2 we denote a set of all subsets of ft, for 
A C ft we denote A' = fl\A. 

A normalized fuzzy set in the set ft is a function A : ft — > [0, 1] such that there exists xq S ft with 
A(xq) = 1. If A(x) < 1 for all x G ft then we call such a function a degenerated fuzzy set. In this 
paper by fuzzy set we always mean a normalized fuzzy set. For any a G (0, 1] an a-cut of a fuzzy set 
A is a set A a = {x G ft : A[x) > a}. 

We will use following notations and definitions. 

• !FS(ft) - Set of all fuzzy sets in a set ft. 

• !FK{ft) - Set of all compact fuzzy sets in a topological space ft, i.e. such fuzzy sets that for 
every a € (0, 1] a set A a is compact and nonempty. 

• J r I(ft) - Set of all fuzzy intervals in a topological and linear space ft. A G TI{ft) if and only if 
A G TJC(ft) and for every a G (0, 1] a set A a is convex. 

• TM{ft) - Set of all fuzzy numbers. A G J 7 J\f(ft) if A is a fuzzy interval and there exists exactly 
one xq G ft such that A(xq) = 1. 

• K,(ft) - Set of all compact subsets of a metric space ft. 

• 2 = /CX(M) - Set of all intervals. 

For x G X we denote ']. Moreover for A G J-T(M.) we define A as a right slope of A, i.e. 

A = A\f^T s. Similarly for a left slope of A. 

Every fuzzy set A in a set ft defines a fuzzy measure II^ on fl defined for any B C ft by 11^4 (B) = 
sup A(B). Any measure obtained in this fashion is called a possibility measure or just a possibility. 
Any such a measure might be defined on a ring of sets T = 2 n containing all subsets of ft, so we may 
not care about its domain. 

Let / : £1 — > E be a function and A G J-S(ft). Then / induces a map /* : TS{£1) — > J-S(E) as 
f*(A)(y) = sup A(f~ 1 (y)). Similarly for any measure 11^4 on ft the function / induces a measure on 
£ as /„(II A )(B) =U A (f- 1 (B)). 
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Let be a set ft with a possibility measure EL A fuzzy variable on f2 with values in a set E 

is any function X : — > £\ Every fuzzy variable induces from II a possibility measure on E called a 
distribution of X. A fuzzy set corresponding to this measure will be called a membership function of 
X. Fuzzy variables having fuzzy intervals as their membership functions will be called fuzzy interval 
variables while those with distributions being fuzzy numbers will be called fuzzy number variables. 

A t-norm is a function T : [0, l] 2 — > [0,1] which is commutative, associative, non-decreasing in 
both variables with 1 being its neutral element. By its properties any i-norm may be extended to 
a function of countable number of arguments. T-norms generalize logical multiplication and serve 
as a definition of independence. We say that B,C C f2 are T-independent for some measure ILa if 
Ha(B P\C) — T(Ua(B), Ha(C)). A generalization to any number of events is routine. Fuzzy variables 
are T-independent if for all B,C C O sets X~ 1 (B) and X~ l {C) are T-independent with obvious 
extension for any number of fuzzy variables. 

2 Convergence 

2.1 Convergence of fuzzy sets 

Let (E, (Ie) be a metric space and denote by 1C(E) a set of all compact subsets of E. The space K,{E) 
is a metric space with a metric d/c(E) defined as 

dic(E)(K,L) = max(sup inf d x (x, y), sup inf d E (x,y)) 

We define a convergence of general compact fuzzy sets as follows. 

Definition 2.1. We say that a sequence (A n )^_ 1 of compact fuzzy sets in a metric space E converges 
to a compact fuzzy set A levelwise if and only if for every a £ (0, 1] we have A" — > A a in K(E). 

In the literature there are used also different types of convergence. Two most important are a 
pointwise convergence of membership functions, e.g. [llj and uniform convergence given by a metric 
d(A, B) = sup Qg ( ji d/c(E)(A a , B a ), e.g. [12] • However, the first one does not agree with philosophy of 
fuzzy sets seeing as a collection of a-cuts, while the latter is too strong if we do not assume compactness 
of supports. Our definition 1 2 . 1 1 will occur to be more natural for fuzzy sets and consistent with fuzzy 
variables convergence. Note, that in this topology maps A i— > A a are continuous for all a S (0, 1]. 
The following proposition links our definition of convergence of compact fuzzy sets with pointwise 
convergence of their membership functions in some special case. 

Proposition 2.1. Let {A n )'^' =1 be a sequence of compact fuzzy sets in a complete metric space E. 

1. If for all a G (0, 1], A" +1 C A" for all n, then (A n )'^' =1 converges to a compact fuzzy set A with 
a-cuts given by A a = H^Li ^n- 

2. If for all a G (0, 1], A" +1 C A" for all n, then for every x G E we have A n (x) — > A(x) with A 
obtained from the previous point. 

3. If for every x G E a sequence (A n (x)) ! ^L 1 is non-increasing and converges to a membership 
function of a compact fuzzy set A, then A n — > A in the sense of the definition \2.1\ 

Proof. 1. By definition we must show that if (K n )^ D =1 is a sequence of compact subsets of E such 
that K n+ i C K n for all n, then (K n )^ =1 converges to K ~ Pl^Li ^"m which is non-empty by 
completeness. Observe that for a given N we have d(K,K^) — sup xeI<N d(K,x). Now assume 
that (K n )^ D =1 does not converge to K, i.e. there exists e > such that for all there exists 
n > N and x G K n with d{x,K) > e. Firstly x £ K, but on the other hand x G K n for all n. 
Contradiction. 

2. By definition A n (x) = sup{a : x G A"} and similarly for A(x). Then, assuming sup0 = we 
have 

oo oo 

A(x) = sup{a : x G C] A"} = sup C] {a : x G A"} = lim sup{a : x G A"} = lim A n (x) 

n— 1 n— 1 

The limit exists since A n (x) is non-increasing. 
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3. It is clear that A^+i C A^ for all a E (0, 1] and hence by point 1 we have A 



n 



A. 



□ 



2.2 Convergence of fuzzy variables 

Firstly, let us recall a definition of fuzzy variable and its membership function. 

Definition 2.2. Let (fl, II) be a set with a possibility measure given by a fuzzy set A in Q. A fuzzy 
variable X with values in a set E is any function X : Q — » E. A possibility measure X*Tl on E 
induced by X is called a distribution of X while a fuzzy set X*A in E determining the distribution is 
called a membership function of X . 

Definition 2.3. Let (X n )^ =1 be a sequence of fuzzy variables defined on a space (fi,II) with values 
in R n (in general in some Polish space). We define following types of convergence. 

1. Almost sure convergence. n({w £ O : linin^oo X n (oj) ^ X(u)}) = 0, 

2. Weak almost sure convergence. II({w £ fl : linin^oo X n (ui) — X(uj)}) = X, 

3. Convergence in measure. For every e > we have liniyj^oo Il(|X n — X\ > e) = 0. 

4- Weak convergence in measure. For every e > we have lim„^oo II(|X n — X\ < e) = 1. 

Moreover if all X n are compact fuzzy variables then we say that [X n )'^L 1 converges to a compact fuzzy 
variable X in distribution if and only if membership functions of X n converge to the a membership 
function of X in the sense of definition 

Observe that a convergence in measure is equivalent to the Fuller's definition [5] when a limit is 
one-point fuzzy variable. We would prefer not to use necessity measures now because we are going 
to discuss them later. Note that a convergence in distribution is similar to a definition given by 
Williamson |27j with pointwise convergence exchanged by levelwise. 

Proposition 2.2. Let (X n )'^' =1 be a sequence of fuzzy variables with values in R™. 

a) There are following dependencies between particular types of convergence: convergence in mea- 
sure =>■ almost sure convergence => weak almost sure convergence => weak convergence in mea- 
sure. 

b) If c is a "non- fuzzy" one-point fuzzy variable, then convergence of distributions is equivalent with 
convergence in measure. 

Proof. For the first arrow we have the following reasoning. Assume for every e > we have liniyi^oo II(|X n — 
X\ > e) = 0. By definition (X„)$£Lj converges almost surely if for any e > we have 



Assume (B n )^ =1 a sequence of subsets of fl such that B n+ \ C B n . Then a sequence (n(B n ))^L 1 is 
non-increasing and has a limit. Moreover II (H^Li Bn) < linin^oo H(B n ). Hence 




< 




N^oo n>N 



lim sup n(|A„ - X\ > e) = 



since by assumption the limit exists. 
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The second arrow is obvious. For the third observe that for every e > we have 

/ OO oo \ 

1 = n u n {\ x ^-x\ 



KN=1 n=N 
/ oo 



□ 



= supn f| {\X n -X\ <e} 

N \Ln J 
< sup lim U(\X n - X\ < e) 

= lim II(|X„ - A| < e) 

n— >oo 

For b) part denote by Ax n membership functions of X n and observe the following 

U(\X n - c\ > e) < 5 4 B CB(c,e) 
where B(c, e) = {x G R" : \x — c\ < e} from which thesis follows. 

3 Tools 

3.1 T-norms, conorms and measures 

The following well known facts can be found in Alsina et al. [2] , Klement [14] or Aczel [1] . 

T-norm T is called Archimedean if for any (0, 1) there exists neN such that T(x, . . . , x) < 

n 

y. Every continuous Archimedean t-norm has a representation T(x,y) = g^^(g(x) + g(y)) where 
g : [0, 1] — > [0, oc] is a strictly decreasing and convex function such that g(l) = called an additive 
generator of T. Here g^~^ : [0, oo] — > [0, 1] is a pseudo-inverse, i.e. g^^{y) = ,9^ 1 (y) for y G [0, g(0)] 
and otherwise. This representation is unique up to a multiplication by a positive constant. A t-norm 
T is called strict if it is strictly increasing. A continuous Archimedean t-norm is strict if and only if 
3 (0) = +cx>. 

If g is an additive generator of a continuous Archemedean t-norm T, than h(x) = e~ g ( x ) is a 
multiplicative generator of T, i.e. T(x,y) — h^~^(h(x)h(y)). Here h : [0,1] — > [0,1] and a psuedo- 
mverse h\~^ : [0,1] -> [0,1] is defined as = /^(y) for y G [h(0),l] and otherwise. A 

continuous t-norm T is Archimedean if and only if it admits a representation by a multiplicative 
generator h : [0,1] — > [0,1]. If /ii and /12 are two generators of T then /i2(x) = and 
ft-i _1 '(y) = /i^ -1 ^") f° r some a > 0. Hence a multiplicative generator is of the form h{x) — e~ a9 ^ 
and so h(l) =1, T is strict if and only if h(0) = 0. 

A t-conorm S, also known as s-norm, is a function 5 : [0, l] 2 — > [0,1] which is commutative, 
associative, non-decreasing in both variables with being its neutral element. By its properties any t- 
conorm may be extended to a function of countable number of arguments. T-conorms generalize logical 
addition and in the theory of decomposable measures serve as a definition of generalized addition. A 
measure fi : T — » [0, 1] is called ^-decomposable and normalized if /i(0) = 0, [i{X) = 1 and for any 
B,C G T such that BCiC = 0, /i(B UC) = S(n(B), A*(C)). A generalization to the countable number 
of sets follows from properties of t-conorms. Further discussion on decomposable measures can be 
found in 18^ and [TS] 

A continuous t-conorm S is called Archimedean if T(x, y) = 1 — 5(1 — x, 1 — y) is an Archimedean 
t-norm. Every continuous Archimedean t-conorm has a representation S(x,y) = ^ (g(x) + g(y)) 
where g : [0, 1] — > [0, 00] is an increasing function such that g(0) = called an additive generator of 
S. This representation is unique up to a multiplication by a positive constant. A t-conorm S is strict 
if it is strictly increasing. A continuous Archimedean t-conorm is strict if and only if g(l) = +00. 

Let A G TS(X) and B G TS(Y). A T-product of fuzzy sets A®B ^ FS{X x Y) for some t-norm 
T is a fuzzy set defined as {A ® B)(x, y) = T(A(x),B(y)). A fuzzy measure coming from A ® B is a 
unique possibility product measure. 
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3.2 Nguyen- Fuller-Keresztfalvi theorem 

If X, Y are two T-independent fuzzy variables with values in E and memberrship functions Ax and Ay 
respectively and / : E x E — > E' is a function, then a fuzzy variable Z = f(X,Y) has a membership 
function given by Zadeh Extension Principle 

A z {z) = sup{T(A x (x), A Y {y)) : f(x, y) = z} 

By Nguyen- Fuller-Keresztfalvi theorem (NFK theorem) [5] if / is continuous, Ax, Ay, T are upper 
semicontinuous and Ax, Ay are compactly supported then 

Az= (J f( A x> A Y) 

with / applied in a set theoretic sense. This theorem can be generalized to the following version 

Theorem 3.1 (Nguyen-Fuller-Keresztfalvi theorem). Let X,Y,Z be topological spaces and let in Z 
every one-point subset be closed. Let f : X xY — > Z be a continuous map, T an upper semicontinuous 
t-norm and A G FK{X), B £ TK{Y). Then 

[MA,B)] a = |J f(AC,B*) 

T(t, v )>a 

for all a G (0, 1] where /» : J-}C(X) x TK,(Y) — > J-)C(Z) is induced from f by Zadeh Extension 
Principle. 

Proof. We need to show three points. 

1. A and B have upper semicontinuous membership functions. Indeed, take a sequence (i n )™ =1 of 
points of X such that x n — > x and (^(xn))^! is non-increasing. Then if for some N G N and 
a G (0, 1] we have xjv G A a , then G A a for all n < N and hence A(x) < limn^oo A(x n ). On 
the other hand by compactness if all x n G A" then a; G A a , so lim n _ >o0 A(x n ) < A(x). 

2. By theorem 1 of [5] we need to show that a function ip(x, y) — {A®B)(x, y) attains its supremum 
on f^ 1 (z) for every z G Z. If this supremum is non-zero then there exists (xo,yo) G f^ 1 (z) such 
that y(a;o, J/o) = e > 0. Because for T(£, 77) > e it is necessary that ^ > e and 77 > e we see that 
the supremum is attained in a set j4 c x B e which is compact. By assumptions A e x B e n / _1 (z) 
is compact and hence y> attains its maximum. 

3. The image of /* lies in TK,{Z). Because by previous point (p attains its supremum, then we 
have by definition 

[f(A,B)] a = {zeZ: sup <p(x,y)>a} = f(<p- 1 (la,l})) 

f(x,y)=z 

A set ( y 9~ 1 ([a, 1]) is a closed subset of a compact set A a x B a and hence [f(A, B)] a is compact. 

□ 

In the formulation of Fuller and Ketereszfalvi the assumption on Z is missing. Now we have a few 
corollaries which are very important in further analysis. 

Corollary 3.2. Let X and Z be topological spaces and let in Z every one-point set be closed. Let 
f : X — * Z be a continuous function and A G J-fC(X). Then for an induced map /* : J-IC(X) — ► J-)C(Z) 
we have {f*A) a = f(A a ). 

Proof. Take for Y a one-point topological space and as B a one-point fuzzy set in Y and apply theorem 

i3~n □ 

Corollary 3.3. I/Ti and T2 are two upper semicontinuous t-norms such that T\ < T2 and X, Y , Z , 
A, B, f as in the NFK theorem\^ then /* (Tl) (A, B) C f { ? 2) {A, B) where f { P stands for a function 
induced from f by means of Zadeh Extension Principle with use of a t-norm T . 
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Proof. If Ti < T 2 then {(£,??) : T x 77) > a} C {(£,77) : T 2 (C,ry) > a} and hence by monotonicity of 
t-norms [/* (Tl) (A, C [/i T2) (A, □ 

Corollary 3.4. Let T be an Archimedean t-norm with an additive generator g and let X , Y , Z , A, 
B, f as in the NFK theorem \3.1\ Then 

[U{A,B)] a = |J 

ff(0+ff(»7)<s(a) 

Proof, g o is not an identity map. If (£, 77) are such that g(£) + gijf) < g(0) then g(£) + g(rj) = 
g (a) ^ T{£, v) = a. If g(Q + g{n) > g(0) then T(f, 77) = but the equation «?(£) + g(r,) < g(a) 
does not have a solution. Hence the only difference is for 0-cut and so is meaningless. □ 

This corollary can be easily extended to a formula valid for any finite number of fuzzy sets. 
3.3 Fuzzy modal value 

A fuzzy modal value is a counterpart of a probabilistic expectation value, however, unlike in probability 
theory, it is a set rather than a number. It is defined as follows. 

Definition 3.1. A fuzzy modal map M : !FS(fl) — > 2 n \{0} C ^5(0) is a function defined as MA = 
A 1 . Similarly for fuzzy variables we define M as the operator with values in 2^\{0} as M.X = M.Ax 
where Ax is a membership function of X . 

If a fuzzy set A is a fuzzy number or a fuzzy variable X is a fuzzy number variable, then we treat 
M.A and MX as a point (number) rather than a one-point set. 

Proposition 3.5. Let f : E — > E' be a continuous function between two topological spaces with E' 
such that all one-point sets are closed. Let A, B € J-K,{E) and X : Q — > E be a fuzzy variable with 
distribution given by A. We can treat MA as a fuzzy set and hence 

1. M{TJC{E)) = K.(E) C TK,(E) 

2. Mf*A = fJAA, 

3. MT(A,B) =T{M(A),M(B)), 

4. Mf(X) = f(MX). 

Proof. By NFK theorem we have the following 

1. by definition, 

2. MfU = {UA) 1 = f{A l ) = /.MA, 

3. MT(A,B) =T(A,B) 1 = (AnB) 1 = A 1 n B 1 =T(A 1 ,B 1 ) = T(MA,MB), 

4. Mf(X) - A) (x) = {UA x f = f{A\) = /(MX). 

□ 

Note that properties 2 and 4 are true for any continuous functions so are more general than their 
counterparts in probability theory because expectation value is a linear operator. 
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3.4 Mean and median 

Let us define two most important functions of type W 1 — > R for any n G N which will br used in the 
further analysis. 

1 ™ 

Ave(xi,...,x n ) = - xj 

U 3=1 

Med(xi, . . . , X2n+i) = median{x\, . . . ,X2n+i) 



Proposition 3.6. Let T be a continuous Archimedean t-norm with an additive generator g. Let 
a G (0,1], n G N and A be a compact fuzzy set in K. Denote N = 2n + 1. TTien 

Med*(A,...,A) Q = ^''(ivTTfM) 

TV 

Proof. The median induces by NFK theorem a map 

Med* (A,..., A)" = |J Med(A^,...,A« N ) 

" ^ ' S(?l) + ---+9(««)<s(a) 

with M ed( A^ 1 , . . . , A^ " ) being such a set among A| which contains and is contained in exactly n 
different sets A^ j . 

Let £i, . . . ,£jv satisfy g(£i) + . . . + .g(£jv) < 3(a)- For any given sets A^ 1 , . . . , A^™ we may order 
them and so assume A^ 1 C . . . C A^ N . Now if we set £x = • • • = £n = 1 an d £n+i ~ £n+i) • • • > Cat = £zv 
then gr(£i) + . . . + g(£ N ) < g(a) and Med(A^, . . . , A ?JY ) = Med(A ? \ . . .,A^ N ). Because for every 
£ < n we have MA C A* 7 C A* , hence Med^ 1 , . . . , A^ N ) is such a set among A^ 1 , . . . , A 4 ™ which is 
contained in all the others. The largest possible such a set requires all A& to be equal with minimal 
£j for j > n + 1 and hence £j- = .g^ 1 ' ^^j<7 (en)) • We obtain thesis by substitution n + 1 = □ 

In the calculation of median we have assumed N to be an odd number. Now we postulate a result 
for any n G N as 

Med* {A, . . . , A) a = A gl ~ 1] (^i^) 



4 Laws of Large Numbers 

Theorem 4.1 (Law of Large Numbers for median). Let T be a continuous Archimedean t-norm with 
an additive generator g and let A G !F]C(M.) For any neN define compact fuzzy sets 

M n = Med* (A, . . . , A) 



Then {M n )^ =l converges to 
Proof. For M„ by proposition 13.61 we have 



MedJA, . . . , A) a = A gl 11 (wttsW) 



JV 



It follows that M£ +1 C M n a and MM n = MA so by propositon[2j]M„ -> M with M Q = f^Li M n 2 
MA. On the other hand if there exists x G K such that x G M" for all n then a; G A^ 3 for all (3 < 1 
and hence a; G MA. Hence (M n )^ x converges to MA. □ 
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Theorem 4.2 (Law of Large Numbers for mean). Let T be a continuous Archimedean t-norm with 
an additive generator g and let A G !FX(M) For any n6N define fuzzy intervals 



A n = Ave* {A,..., A) 



Then (A„)£° =1 converges to MA. 



Proof. First assume A is continuous and concave, i.e. A and A are continuous and concave functions. 
Then by means of works [?], 0, [H], [H>] or [22] we find 

Ave*{A, . . . , A) a =A^ 1] ^ 3{a ^) 



Assume the right slope A of A is not concave and a = A 1 . Denote 0(a) = A a and note that 
: (0, 1] — > [a, oo) is a left-continuous non-increasing function with 0(1) — a. Now set a £ (0, 1] and 
denote Y = (0, 1] x [a, oo) and 

C = {(fl,x) S Y : x < 4>{a)} 

and C — conv (C) a convex hull of C. For any Xq G [a, oo) denote C XQ — {/3 G (0, 1] : (/3, Xq) G C}. 
Note that C is closed in Y, so is C and C^o ^r any xq G [a, oo) and hence denote H(x) — maxC^. 
Moreover H{x) > H(y) if x < y and by definition if (1, x) G C then x — a and hence for any a; G (a, oo) 
we have 7?(x) < 1. H is continuous since any line passing through {x,H(x)) and (<p(a),a) meets C y 
in some point and 

\H(x) - H(y) = maxC^ - max6 a < — — r-\x - y\ 

cp(a) — max(s, y) 

Now we can define a fuzzy interval H : [a, oo) — ► [0, 1] as follows 



H(x) if x < 0(a) 
t4(x) if a; > 0(a) 



_ff defines a fuzzy interval on [a, oo) with continuous membership function which is concave, and in 
particular logarithmically concave on [a, 0(a)]. 

Now carry out the same procedure for the second slope of A, concatenate the fuzzy sets in order 
to obtain a fuzzy interval B on K with continuous membership function and concave on A a . Observe 
that A C B in order to get for any a G (0, 1] in a similar way as in the proof of theorem 14. II 

K C B% 1/n MB = MA 
and finish the proof. □ 



By proposition 12.21 for the convergence of fuzzy number variables it is enough to prove the conver- 
gence of distributions. Hence we have 

Theorem 4.3 (Strong Law of Large Numbers). Let T be a continuous Archimedean t-norm and let 
A G J-J\(M). Let (X n )^ =1 be a sequence of T -independent fuzzy number variables defined on a space 
(0,11) with distributions given by A. For any n G N define fuzzy number variables 

A n = Ave{X ll ...,X n ) 
M n = Med(X l7 ...,X n ) 

Then and (M^^L-^ converge in distribution, in measure, and almost surely to MA viewed as 

a one-point fuzzy variable. 

This theorem was proved by Triesch in paper [23] . However, it was neither interpreted in a context 
of fuzzy variables nor compared with dependencies of fuzzy sets. 
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5 Set complements, necessity measures and logics 



Consider the following facts. 

• We may treat fuzzy sets as a dual theory to probability by means of the formula P(A) = 
J A(x)dP(x) where A is a fuzzy set and P a probability measure. However, a definition of 
Lebesgue-like integral involves addition and multiplication, so a choice of a t-norm should define 
a choice of a f-conorm representing addition. 

• For a given possibility measure II there exists a necessity measure II' given by a formula li'(A) — 
1 — 11(^4'). However this definition should depend on a t-norm since for example a relation 
A' (~l B = B\A involves t-norm. 

• In the sense of necessity it is reasonable to say that a possibility of not A is equal to 1 — H{A) 
and probably not identical with H(A'). 

5.1 Logics and norms 

We would like to define a logic-like structure of a type ([0, 1], ©, <g>, -i) where ffi, © and -i stand for 
logical or, and and negation. Moreover we would like to say that a measure /i has its values in this 
structure in such a sense that 

1. if A n B = then fi(A U B) = © At(-B), i.e., /i is ©-decomposable measure, 

2. if n(A fl B) = fJ,(A) © n(B) then we say that A and B are independent, 

3. n{A') = -i/z(j4), where A' = X\A. 

Due to results of 2 it is impossible for this structure to be a Boolean logic, since if © and © are 
mutually distributive then © = max and © = min. However we can tackle with this problem partially 
if we observe that it is sufficient to define © on a domain D C [0, l] 2 such that (x, y) £ D if and only 
if x © y < 1 or, by theorems of Aczel jT], extend S and T to [0, oo). 

On the other hand observe that in such a logic laws involving sums should be satisfied only if their 
arguments come from measures of disjoint sets. For example de Morgan laws are too general since 
they describe a situation of arbitrary sets. For example in a very natural structure ([0, 1], + , •, 1 — id) 
which leads to the probability theory de Morgan laws are not satisfied. 

Looking for a condition involving © which comes only from disjoint sets we find a following one 

(AnB)U(AnB') = A 

which leads to the definitions 

Definition 5.1. A strict negation n is a strictly decreasing function n : [0, 1] — -> [0, 1] such that 
n(0) = 1 and n(l) = 0. 

Definition 5.2. A normal triple (S,T,n) is a triple consisting of a t-conorm S, t-norm T and a 
strict negation n such that 

S(T(x,y),T(x,n(y))) = x 

for all x,y £ [0, l] 2 . We say that a normal triple (S,T,n) is continuous if S, T and n are continuous 
functions. 

Theorem 5.1. If (S,T,n) is a continuous normal triple, then there exists a unique continuous and 
strictly increasing function h : [0, 1] — > [0, 1] such that h(0) = 0, h(l) = 1 and 

S(x,y) = h^(h(x)+h(y)) 
T{x,y) = h-\h{x)h{y)) 
n{x) = br\l-h{x)) 

In particular T is a strict continuous Archimedean t-norm and S is a strict continuous Archimedean 
t-conorm. 

Proof. This is theorem 3.2.13 in □ 
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5.2 Measures and fuzzy sets 

Now wc postulate the following. Define two fields (cr-fields) T and T' of subsets of such that if 
A E T then A' G T' and vice versa. Sets contained in T will be called events, while sets in T' 
anti-events. 

A function A : (f2, J-h, J-jfj) — > (E,Te,F'e) is called measurable or a fuzzy variable if for every 
A G Te and i? G ^ we have A" 1 ^) G and X" 1 ^) G T[ v 

For a given strict continuous Archimedean i-norm T we take a unique i-conorm 5 and a strict 
negation n such that (S, T, n) constitues a continuous normal triple. For a possibility measure IT : 
T — > [0, 1] we define a necessity-like measure II' : T 1 — ► [0, 1] by II'(A') = n(n(^4)). Then a following 
diagram commutes 

(^, 0, A, C, U, n) ([o, l], 0, 1, <, max, T) 

(-)' 

, 0, A, c, u, n) ([o, l], o, 1, <, S, T) 

We interpret it as follows. If we know a possibility measure of some event A than we can define a 
measure of not A by W(A'). However if A' G T then typically IL(A') ± W(A'). 

For A G T and B E T' sets iUB and An B might not belong neither to T nor T' . However, if 
AnB = 0we can define U(A UB) = max(n(A), II' (2?)) and similarly if A and £?' arc T-independent, 
then U(A DB)= T (11(A), U'(B)). 

Finally, some membership functions of fuzzy variables might become degenerated i.e. sup A(E) 
may be less than 1 and should be normalized. 

Definition 5.3. Let T be a strict continuous Archimedean t-norm with a multiplicative generator h. 
Let A be a degenerated or not fuzzy set in E. A normalized version of A denoted by A N is a fuzzy set 
defined as 

In general m&x(A N , B N ) ^ max(A, B) N . All fuzzy sets should be normalized. 

Here we have to mention that this procedure is not well defined and ambiguous. We will show how 
it works in the proof of Glivenko-Cantelli theorem but we do not know how to tackle with problem 
of normalization thoroughly. 

6 Theory of fuzzy measurement 
6.1 Combination of errors 

Proposition 6.1. Let T be a continuous Archimedean t-norm and let A G FAf(M). Let (A„)^L 1 be a 
sequence of T -independent fuzzy number variables all with distributions given by A. Let Y be a fuzzy 
variable with a distribution given by a characteristic function of some set B and arbitrary dependence 
with X n and let Z n = X n + Y for all n G N. Define 

A n = Ave(Z 1 ,...,Z n ) 
M n = Med(Z u ...,Z n ) 

Then (A n )^ =1 and (M„)^L 1 converge in measure and almost surely to Y. 
Proof. Observe that for any x\, . . . , x n , y G R there is 

Ave(xi+y,...,x n +y) = Ave(xi, . . . ,x n ) + y 
Med(xi+y,...,x n +y) = Med(x\, . . . , x n ) + y 
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which is also true for intervals. Hence by NFK formula A n = A n + Y and M n = M n + Y where 

An = Ave(X 1 ,... 1 X n ) 
M n = Med{Xx,...,Xn) 

By theorem 14.31 we know that A n and M n converge to MA in measure. By 

TLQAn - A\ > e) = n(|l„ - A\ > e) n -^° 

we obtain A n — > A in measure. The same holds for M n . □ 

This proposition allows to treat both random and systematic components of uncertainty in a 
consistent way in the fuzzy variables theory. A fuzzy number represents a random component and is 
obtained due to carrying out a measurement. A non- fuzzy set represents a systematic component of 
uncertainty and is obtained by use of expert method. Such an approach with an arithmetic based on 
i-norms was proposed in [24) and [25) and studied empirically in [26] . 

The following section develops a theorem necessary for a construction of empirical membership 
function of a fuzzy process. 



6.2 Estimation of membership function 



Now we will prove a fuzzy version of Glivenko-Cantelli theorem. Assume we are given a strict 
continuous Archimedean i-norm T with a multiplicative generator h. By theorem 15.11 there exists 
unique i-conorm S and a strict negation n such that (T, S 1 , n) is a normal continuous triple. 

For E = W 1 we consider a held of events T containing all closed subsets of R™ . Similarly a field of 
anti-events T' contains all open subsets of R™. For example if we are given two T- independent fuzzy 
variables X, Y : — > M and two points then the expression II(A < x, Y < y) is equal 

U(X <x,Y <y) = T(U(X < x), U'{Y < y)) = T(U(X < x),n{U(Y > y))) 

If Xi, . . . ,x n is a sequence of real numbers such that x\ < xi < . . . < x n , then denote Xk :n = £fc- 
If x%, . . . , x n is any sequence of real numbers, then Xk :n is equal Xk after ordering. Obviously x\. n — 
min(a;i, . . . , x n ) and x n:n = max(xi, . . . , x n ). Under these assumptions and notations we have 

Theorem 6.2 (Glivenko-Cantelli theorem). Let T be a strict continuous Archimedean t-norm with 
a multiplicative generator h and (X n )^_ 1 be a sequence of T -independent identically distributed fuzzy 
number variables with a continuous membership function A G J-J\f(M). Denote <p(a) = A a and define 
estimators of (j) as 



4>n{a) = 



X rn h ( a ) 



V 



,x 



-I 



j:r, 



if n is even 
if n is odd 



Then 

1. M.cj) n (a) n ^ff (f}(a) for every a 6 (0, 1], 

2. for every e > 

D n = sup d(4>n(a),<t>(a)) 

o£[e,l] 

converges to zero in measure and almost surely. 

The idea of building a membership function of the process is presented in the figure [T] Unfortu- 
nately, the estimators <p n are slightly biased. In order to analyze this bias we have to prove a more 
technical version of the theorem. 

Theorem 6.3 (Glivenko-Cantelli theorem). Let T be a strict Archimedean t-norm with a multiplica- 
tive generator h and (Xn)n°=i be a sequence of T -independent identically distributed fuzzy number 
variables with a continuous membership function A G J-J\f(M.). Let b : N x (0, 1] — ► K be any function 
such that 
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Figure 1: An example of a fuzzy histogram obtained from 8 measurement points. 

1. limn^oo b(n, a) — for all a G (0, 1], 

2. k{n,a) = |(/i(a) + b(n, a)) € {1,2,..., [§]} /or a!/ n G N anda G (0,1], 
5. if a < (3 then k{n 1 a) < fc(n, /3) /or a^/ n G N. 

Finally denote 4>{a) — A a and define estimators of (f> 

0n(oO = \Xk(n,a)\m -^-n—k(n ,a):n] 

Then 

1. Mcj) n (a) = (t>(hr l (h(a) + b(n,a))) for every a G (0,1], 

2. for every e > 

£>„ = sup d((f) n {a),<j){a)) 

a£[e,l] 

converges to zero in measure and almost surely. 

In this formulation b is a bias which can be added by hand and we see that <f> n are quite immune 
to bias. From this theorem we have 

Corollary 6.4. Theorem 1 6. 2\ holds with \b(n,a)\ < 2/n. 

Proof of the theorem \6.'A For a first point set n G N and a G (0,1] and take any x — [x,x] G T. 
Denote k = k{n, a) and observe that 

H(4> n (a) = x) = n(there are exactly k indices such that Xj < x and k such that Xj > x) 



) > x, 



II I [J {X a (i) < X,. . . , X^k) < ^a(n-k+l) > X,.. ., X a ( n 

{X a ( k +i) <x\J X a{k+1 ) > x)', . . . , {X a{n _ k) <iV X ff ( n _ fe ) > x)'}) 
n (Xi < x, . . . , X k < x, X n - k+ i >x,...,X n >x 

(X k+1 <xV X k+1 > x)', (I n - fc <iV X n _ k > x)') 

= t I u(Xi <x),...,n(Xi <x),n(X! >x),...,n(Xi 

k k 



n-(n(Xi <iVli >#)),..., n(n(Xi <iVli > x)) 

n-2fc 

h,- 1 (h k (A(x)) ■ h k {A{x)) • (1 - h{mzx{A{x),A{x)))) n - 2k 
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Since h is an increasing homeomorphism we may set y — h(A(x)) and z — h(A(x)) and proceed as 
follows. A degenerated fuzzy set B n (y, z) — (y k z k (l — max(y, z)) n ~ 2k ) defined on [0, l] 2 has a unique 
maximum at yo = zq = — for y, z € [0, l] 2 and by continuity of A and monotonicity of h there exists 
unique x = h (4 > (2k/n)) at which maximum is attained. Hence 

M<f> n (a) = 4> (h- l {h{a) + b(n, a))) 

For a second part by assumption 1 on 6 and continuity of h and <j> we obtain M<j> n {a) — + A a for 
all a€ (0,1]. A normalization of B n gives 

(j, z )(M«)+6(».«))/2(i - max(y, z ))i-H<x)-b(n,a) 
{h{a) + b{n, a ))H<*)+b(n, a ) J{ _ h ^ _ b ^ a ))HW-H»:«) 

By symmetry we may limit ourselves to the case z < y. Substituting y n = h(a) + b(n, a) we find 
{yn-, Vn) = 1- Increasing B^ by taking y — z we obtain 



B l : < F n (y) 



yh(a)+b(n,a) M _ y\ 1— h(a)— b(n,a) 



,(/i(a) + b{n,a)) h ( a )+ b ( n > a )(l - h{a) - b(n, a)y- h (»)-b(n, a ) 
Taking a derivative of F n we have 

dhxF n y — h(a) — b(n,a) 
dy y{y - 1) 

so F n < 1 for all y G [0, l]\{y n } which means that B^ ^—3? l^h(a).h(a))} pointwise for all a S (0, 1]. 
For monotonous convergence take B™ (Hi z ) = su Pfc>n &k iVi z ) an d observe that fuzzy sets B^ are 
subsets of B!% . Since B^ (Vi z ) ^S 3 lr(h(a)M a ))} P om twise, then by a proposition ^. II <f> n (a) — > 0(a) 
in distribution. 

For uniform convergence apply the standard procedure. We will focus on the left slope because 
the procedure for the right one is the same. Take e > and Mgl Define points 

Ae A e -A\ 
XM,k = A H — — k 

for k = 0, 1, . . . , M. and ak = A(xk). If a is such that a £ [ak-, ctfc+i] then we have 

±Jp) ~ 0(a) < n ( a fe+i) ~ #*fc) < ^("fe+i) " + 
as well as <j> (a) — 4>(a) > 4> n (ak) — 4>{a-k) — -jj- Hence we have 

\<t n {a) -0(a)] < max(0 ji (a fc+ i) - 0(a fc+ i),0 n (a fc ) - 0(a fc )) + i = max(A^^ 
where we defined two As for simplicity. Thus we have 

1} = max( max Ai + h n max A1°m *•) + T7 



A (0) ) + — 

\=n,M,k 1 =n,M,k> 1 ^ 



We have shown that any A converges to zero, hence finite maxima converge to zero. Thus we have 
limsuPrc ,^ D n < 1/M. The same reasoning for D n completes a proof. 

□ 

Note that if F : M — > [0, 1] is a continuous cumulative distribution of a random variable A then a 
corresponding fuzzy distribution A 6 JFA/"(]R) is given by 



A" = 



supF- A ( -J, infill- 



This is a probability-possibility transformation similar to those presented in [3] such that A a = I a 
where I a is a a-confidence interval defined as the right hand side of the equation above. 
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6.3 Realization of fuzzy variables 



In the theory of probability we can simulate realizations of random variables by means of inverse 
distribution function. A similar theorem holds for fuzzy variables. 

Proposition 6.5. Let U : Q — > [0,2] be a fuzzy variable with a membership function Au(x) — 
a;l[o + (2 — 01)1(1^] (i). Let X be a fuzzy number variable with a continuous membership function 
A x e FN(B) an & define a function t/> : [0, 2] — ► K 



ip(ai) 



A' 1 (a) if a < 1 



Ax (2 -a) if a>l 
This function is continuous and ifi(U) = X almost surely. 

Proof. We have TL(ip(U) = x) = U.(U = il)~ l {x)) which for x < MA X is equal U(U = A(x)) = A(x). 
For x > MA X we have Il(U = ip~ 1 (x)) = IL(U = 2- A(x)) = A(x). □ 

This proposition is a generalization of probability-possibility transformations and shows that only 
fuzzy numbers represent probabilistic quantities. It allows to simulate fuzzy processes and, together 
with Glivenko-Cantelli theorem and the Strong Law of Large Numbers shows, that the theory is 
consistent in this sense, that it allows to produce realizations of fuzzy variables and then regain their 
distributions by means of membership function estimators. 



7 Conclusions and generalizations 

Theorem 16.31 allows to construct a membership function of a fuzzy process from empirical data by a 
direct estimation of its a-cuts. This theorem guarantees that the presented procedure described in 
the theorem converges to a real membership function. 

If empirical data are clearly probabilistic, then a multiplication is a good choice of a i-norm. How- 
ever for real processes another choice of a i-norm might be more suitable. The problem of estimation 
of a t-norm was raised and analyzed empirically in |26j . The procedure is based on minimalization 
of a distance between fuzzy sets obtained by empirical averaging and theoretical one with use of a 
i-norm. A distance of two fuzzy sets is given by 

d(A,B) = sup d K{x) {A a ,B a ) 

t»e[e,l] 

for some set e > 0. 

Now we must point out some important problems encountered in the procedure of normalization 
and proof of Glivenko-Cantelli theorem. The procedure of normalization developed in this paper 
seems to be ambiguous. The second problem is encountered when a generalization of Glivenko- 
Cantelli theorem is considered. The assumption about continuity of a membership function of a 
process is important and cannot be rejected. In order to make a possibility of use non-strict i-norms 
we should weaken a definition of strict negation to the following one. 

Definition 7.1. A negation is a function n : [0,1] — > [0,1] which is continuous, strictly decreasing 
and such that n(l) = 0. 

With this set-up it is possible to generalize theorem 15.11 and obtain all Archimedean i-norms. 
However then the proof of Glivenko-Cantelli theorem breaks down since h is not a homeomorphism. 
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